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Abstract 

Electronic transport characteristics through a single phenalenyl molecule sandwiched between two metallic 
electrodes are investigated by the use of Green's function technique. A parametric approach, based on 
the tight-binding model, is used to study the transport characteristics through such molecular bridge 
system. The electronic transport properties are significantly influenced by (a) the molecule-to-electrode 
interface structure and (b) the molecule-to-electrode coupling strength. 
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1 Introduction 

Molecular transport have attracted much more at- 
tention since molecules constitute promising build- 
ing blocks for future generation of electronic de- 
vices. Electron transport through molecules was 
first studied theoretically in 1974 by Aviram et 
al. [T]. Since then numerous experiments [2l [3l 
HI m [6] have been performed through molecules 
placed between two electrodes with few nanome- 
ter separation. The operation of such two-terminal 
devices is due to an applied bias. Current pass- 
ing across the junction is strongly nonlinear func- 
tion of the applied bias voltage and its detailed 
description is a very complex problem. The com- 
plete knowledge of the conduction mechanism in 
this scale is not well understood even today. The 
transport properties of these systems are associ- 
ated with some quantum effects, like as quantiza- 
tion of energy levels, quantum interference of elec- 
tron waves [a[8l[9l[T0l[ni[12l[13l[l4],etc. Follow- 
ing experimental developments, theory can play a 
major role in understanding the new mechanisms 
of conductance. The single-molecule electronics 
plays a key role in the design of future nanoelec- 
tronic circuits, but, the goal of developing a re- 
liable molecular-electronics technology is still over 
the horizon and many key problems, such as device 
stability, reproducibility and the control of single- 
molecule transport need to be solved. It is very es- 
sential to control electron conduction through such 
quantum devices and the present understanding 
about it is quite limited. For example, it is not 
very clear how the molecular transport is affected 
by the structure of the molecule itself or by the na- 
ture of its coupling to the electrodes. To design 
molecular electronic devices with specific proper- 
ties, structure-conductance relationships are needed 
and in a recent work Ernzerhof et al. [15j have pre- 
sented a general design principle and performed sev- 
eral model calculations to demonstrate the concept. 

There exist several ah initio methods for the cal- 
culation of conductance [Ml [HI [HI [HI EHl Ell [22] 
through a molecular bridge system. At the same 
time the tight-binding model has been extensively 
studied in the literature and it has also been ex- 
tended to DFT transport calculations [23]. The 
study of static density functional theory (DFT) [21] 
within the local-density approximation (LDA) to in- 
vestigate the electronic transport through nanoscale 
conductors, like atomic-scale point contacts, has 
met with nice success. But, when this similar theory 
applies to molecular junctions, theoretical conduc- 



tances achieve larger values compared to the ex- 
perimental predictions and these quantitative dis- 
crepancies need extensive study in this particular 
field. In a recent work, Sai et al. [25] have predicted 
a correction to the conductance using the time- 
dependent current-density functional theory since 
the dynamical effects give significant contribution in 
the electron transport, and illustrated some impor- 
tant results with specific examples. Similar dynam- 
ical effects have also been reported in some other re- 
cent papers [261 El] , where authors have abandoned 
the infinite reservoirs, as originally introduced by 
Landauer, and considered two large but finite op- 
positely charged electrodes connected by a nano- 
j unction. 

Our aim of the present article is to reproduce an 
analytic approach based on the tight-binding model 
to characterize the electronic transport properties 
for the model of a single phenalenyl molecule and fo- 
cus our attention on the effects of (a) the molecule- 
to-electrode coupling strength and (b) the quantum 
interferences in these transport. This is an interest- 
ing molecular system to study because of the unique 
behavior of the phenalenyl molecule and contact- 
ing the molecule to the electrodes from two differ- 
ent locations several interesting results are obtained 
for the electron transport. Here we adopt a simple 
parametric approach [281 ES [30l [311 [32l [33] for this 
calculation. The parametric study is motivated by 
the fact that the ab initio theories are computation- 
ally very expensive and here we concentrate only on 
the qualitative effects rather than the quantitative 
ones. This is why we restrict our calculations only 
on the simple analytical formulation of the trans- 
port problem. Not only that, the ab initio theories 
do not give any new qualitative behavior for this 
particular study in which we concentrate ourselves. 

This paper is specifically arranged as follows. In 
Section 2, we give a very brief description for the 
calculation of transmission probability (T), con- 
ductance {g) and current (/) through a finite size 
conductor sandwiched between two metallic elec- 
trodes. Section 3 focuses the behavior of the 
conductance-energy and the current-voltage char- 
acteristics for the single phenalenyl molecule and 
studies the effects of (a) the molecule-to-electrode 
interface structure and (b) the molecular coupling 
strength on the above mentioned characteristics. 
These two factors i.e., the interface structure and 
the coupling strength play significant role on the 
quantum transport through molecular devices. Fi- 
nally, we summarize our results in Section 4. 
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2 A glimpse onto the theoret- 
ical formulation 

Here we describe very briefly about the methodol- 
ogy for the calculation of transmission probability 
(T), conductance (g) and current (/) through a fi- 
nite size conducting system attached to two semi- 
infinite metallic electrodes by using the Green's 
function technique. 

Let us first consider a one-dimensional conduc- 
tor with number of atomic sites (array of filled 
circles) connected to two semi-infinite metallic elec- 
trodes, namely, source and drain, as given in Fig.[T] 
The conducting system in between the two elec- 
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Figure 1 : Schematic view of a one-dimensional con- 
ductor with N number of atomic sites (filled cir- 
cles) attached to two electrodes, source and drain, 
through sites 1 and A^, respectively. 

trodes can be an array of few quantum dots, or a 
single molecule, or an array of few molecules, etc. 
At low voltages and temperatures, the conductance 
of the conductor can be written by using the Lan- 
dauer conductance formula [34l [35] , 
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where g is the conductance and T is the transmis- 
sion probability of an electron through the conduc- 
tor. The transmission probability can be expressed 
in terms of the Green's function of the conductor 
and the coupling of the conductor to the two elec- 
trodes by the expression [SU [35] , 

(2) 



T = Tr [TsG^cTdG^c] 



where and are respectively the retarded and 
advanced Green's function of the conductor. Fg 
and To are the coupling terms due to the coupling 
of the conductor to the source and drain, respec- 
tively. For the complete system i.e., the conductor 
including the two electrodes, the Green's function 
is defined as. 



G={e-H) 



(3) 



where e = E -\- irj. E is the injecting energy of the 
source electron and rj gives an infinitesimal imagi- 
nary part to e. Evaluation of this Green's function 



requires the inversion of an infinite matrix as the 
system consists of the finite conductor and the two 
semi-infinite electrodes. However, the entire system 
can be partitioned into sub-matrices corresponding 
to the individual sub-systems, and the Green's func- 
tion for the conductor can be effectively written as. 



Gc — (e — He — — 



(4) 



where He is the Hamiltonian for the conductor 
sandwiched between the two electrodes. The single 
band tight-binding Hamiltonian for the conductor 
within the non-interacting picture can be written in 
the form, 



He = ^iclci 



<ij> 



(5) 



where c| (c^) is the creation (annihilation) opera- 
tor of an electron at site i, e^'s are the site energies 
and t is the nearest-neighbor hopping integral. In 
Eq. [H S5 = h\(jgshsc and S^, hDcgoh^Dc are 
the self-energy operators due to the two electrodes, 
where gs and go are respectively the Green's func- 
tion for the source and the drain, hsc and hjjc are 
the coupling matrices and they will be non-zero only 
for the adjacent points in the conductor, 1 and N 
as shown in Fig. [1] and the electrodes respectively. 
The coupling terms Fg and F^i for the conductor 
can be calculated through the expression. 
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where E^^ ^-j and ^-j are the retarded and ad- 
vanced self-energies, respectively, and they arc con- 
jugate with each other. Datta et al. [34] have shown 
that the self-energies can be expressed like as, 
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(7) 



where ^{s.d} are the real parts of the self-energies 
which correspond to the shift of the energy eigen- 
values of the conductor and the imaginary parts 
^{S,D} of the self-energies represent the broaden- 
ing of these energy levels. This broadening is much 
larger than the thermal broadening and this is why 
we restrict our all calculations in this article only at 
absolute zero temperature. The real and imaginary 
parts of the self-energies can be determined in terms 
of the hopping integral {t^s,d}) between the bound- 
ary site of the conductor and the electrodes, the in- 
jection energy (E) of the transmitting electron and 
the hopping strength (v) between nearest-neighbor 
sites of the electrodes. 
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Thus the couphng terms Ts and To can be writ- 
ten in terms of the retarded self-energy as, 



{S,D} 



= -2Im 



^{S.D} 



(8) 



Now all the information regarding the conductor 
to electrodes coupling are included into the two 
self energies as stated above and are analyzed 
through the use of Ncwns- Anderson chemisorption 
theory [551 122] ■ The detailed description of this the- 
ory is obtained in these two references. 

Hence, by calculating the self-energies, the cou- 
pling terms Ts and F^i can be easily obtained and 
then the transmission probability (T) will be com- 
puted from the expression as mentioned in Eq. [21 

Since the coupling matrices hsc said hoc sltc 
non-zero only for the adjacent points in the conduc- 
tor, 1 and N as shown in Fig. [1] the transmission 
probability becomes, 

T{E, V) = 4Afi(i;, V)A^^iE, V)\GiNiE, V)\^ 

(9) 

=< iV|A|iV > and 



where An =< 1|A|1 >, 
GiN^< l\Gc\N>. 

The current passing through the conductor is de- 
picted as a single-electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current-voltage relation is evaluated from the fol- 
lowing expression [35] . 



eV/2 
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T{E, V)dE 



(10) 



EF-eV/2 



where Ep is the equilibrium Fermi energy. For the 
sake of simplicity, here we assume that the entire 
voltage is dropped across the conductor-electrode 
interfaces and this assumption does not greatly 
change the qualitative behaviors of the I-V char- 
acteristics. This assumption is based on the fact 
that the electric field inside the molecule, especially 
for short molecules, seems to have a minimal effect 
on the conductance-voltage characteristics. On the 
other hand for quite longer molecules and high bias 
voltage, the electric field inside the molecule may 
play a more significant role depending on the inter- 
nal structure of the molecule [34|, yet the effect is 
too small. Using the expression of T{E, V) as in 
Eq. [ithe final form of I{V) will be, 

EF+eV/2 

m = ^ J Af,iE,V)A^^iE,V) 

Ep-eV/2 

x\GiN{E,V)\^dE (11) 



Eqs. m [21 and [TTl are the final working formulc for 
the calculation of conductance g, transmission prob- 
ability T, and current /, respectively through any 
finite size conductor sandwiched between two elec- 
trodes. 

By using the above methodology, in this article 
we will investigate the electronic transport char- 
acteristics through a single phenalenyl molecule 
(Fig. 12). Throughout this article we set Ep = 
and choose the unit c — e ~ h ^ 1. 

3 Results and their interpre- 
tation 

In this section wc focus on the conductance- 
energy and current-voltage characteristics of a sin- 
gle phenalenyl molecule and investigate the depen- 
dence of these characteristics on (a) the molecule- 
to-electrode interface structure and (b) the molec- 
ular coupling strength. The schematic representa- 
tions of single phenalenyl molecules attached to the 
two metallic electrodes are shown in Fig. O To 
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Figure 2: Schematic view of single phenalenyl 
molecules attached to two metallic electrodes, 
namely, source and drain, through thiol (S-H) 
groups. 

characterize the molecule-to-electrode interface ef- 
fect, here we consider two different geometries (see 
Figs.j^a) and (b)). These single molecules are con- 
nected to the electrodes by thiol (S-H) groups. In 
experiments, two electrodes made from gold (Au) 
are used and molecules attached to the electrodes 
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by thiol (S-H) groups in the chcmisorption tech- 
nique where hydrogen (H) atoms remove and sulfur 
(S) atoms reside. Here the molecule is described by 
the similar tight-binding Hamiltonian as prescribed 
in Eq. [5] Throughout the article we describe all the 
essential features of electron transport in two dis- 
tinct regimes. One is t^s,d} «t, called the weak- 
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Figure 3: Conductance 3 as a function of the en- 
ergy E for the phenalenyl molecule connected to 
the electrodes according to Fig. [2l^a). (a) and (b) 
correspond to the results for the weak and strong 
molecule-to-electrode coupling limits, respectively. 

coupling limit and the other one is t^s,d} t, called 
the strong-coupling limit, where ts and td are the 
hopping strengths of the molecule to the source and 
drain, respectively. The common set of values of the 
parameters used in our calculations for these two 
limiting cases arc: ts = tjj = 0.5, t = 2.5 (weak- 
coupling) and Ts — Tjj ^ 1.5, t ~ 2.5 (strong- 
coupling). We set the nearest-neighbor hopping 
strength v = 4. 

Figure [3] shows the variation of the conductance g 
as a function of the injecting electron energy E for 
the molecular bridge given in Fig. [2l^a). Figure [S]^ a) 
corresponds to the result for the bridge system in 
the weak coupling limit. The conductance is almost 
everywhere zero, except at the resonances where it 
approaches to 2. The resonant peaks in the con- 
ductance spectrum coincide with the energy eigen- 



values of the single phenalenyl molecule. Therefore 
the conductance spectrum manifests itself the elec- 
tronic structure of the molecule. On the other hand, 
in the strong coupling limit the resonant peaks get 
substantial widths, as shown in Fig.[3l^b) and it em- 
phasizes that the electron conduction takes place 
throughout the energy range (the range taken here 
in the figure). This is due to the broadening of the 
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Figure 4: Conductance g as a function of the en- 
ergy E for the phenalenyl molecule connected to 
the electrodes according to Fig. [^^b). (a) and (b) 
correspond to the results for the weak and strong 
molccule-to-clcctrode coupling limits, respectively. 

molecular energy levels in the strong coupling case, 
where the contribution comes from the imaginary 
parts of the self-energies as described earlier. 

To characterize the molecule-to-electrode inter- 
face effect on electron transport, we plot the con- 
ductance in Fig.[¥]for the molecular bridge given in 
Fig. ^b). Figures SKa) and (b) correspond to the 
results for the bridge system in the weak and strong 
coupling limits, respectively. The broadening of the 
conductance peaks in the strong coupling limit (see 
Fig. H^b)) is due to the same reason as mentioned 
earlier. From the curves plotted in Fig. IH we see 
that the conductance peaks do not reach to unity 
anymore and get much reduced value, compared to 
the results described in Fig. [31 Such a behavior can 
be understood as follows. The electrons are car- 
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ricd from the source to drain through the molecule 
and thus the electron waves propagating along the 
two arms of the molecular ring may suffer a phase 
shift between themselves, according to the result of 
quantum interference between the various pathways 
that the electron can take. Therefore, the probabil- 
ity amplitude of the electron across the molecular 
ring becomes strengthened or weakened (from the 
standard interpretation of the wave functions). It 
emphasizes itself especially as transmittance can- 
cellations and anti-resonances in the transmission 
(conductance) spectrum. Here the mentioned phase 
shift is observed by the variation of the geometry of 
the molecular bridge. The most significant issue for 
considering these two different molecular bridge sys- 
tems is that in this way the interference conditions 
are changed quite significantly. 

Another key feature observed from the conduc- 
tance spectrum is the existence of the conductance 
(transmittance) zero. From Fig. lU it is observed 
that the conductance zero appears across E = 0. 
Such anti-resonant state is specific to the interfero- 
metric nature of the scattering and does not occur 
in usual one-dimensional scattering problems of po- 
tential barriers. It is also observed that the posi- 
tion of the anti-resonant state on the energy scale 
is independent of the molecule-to-electrode coupling 
strength. The width of this state is very small and 
hence it does not give any significant contribution 
to the current-voltage characteristics. However, the 
variations of the interference conditions have strong 
influence on the magnitude of the current flowing 
through the bridge systems. 

Thus it can be emphasized that the electron 
transmission is strongly influenced by the molecule- 
to-electrode interface structure. 

The scenario of electron transfer through the 
molecular junction is much more clearly visible 
from the current-voltage characteristics. Current 
through the molecular system is computed by the 
integration procedure (given in Eq. [11]) [35] of the 
transmission function T. The behavior of the trans- 
mission function is similar to that of the conduc- 
tance variation since g = 2T (from the Landauer 
conductance formula, Eq. [T] with e = h ~ 1 
in our present formulation). In Fig. [5] we plot 
the current-voltage characteristics for the molecu- 
lar bridge shown in Fig. [^J^a). Figures [Sl^a) and (b) 
correspond to the currents for the bridge system in 
the weak and strong coupling cases, respectively. 
It is observed that in the weak coupling case the 
current shows staircase-like structure with sharp 
steps. This is due to the discreteness of molecu- 



lar resonances as shown in Fig. ^a,) . As the volt- 
age increases, the electrochemical potentials on the 
electrodes are shifted and eventually cross one of 
the molecular energy levels. Accordingly, a cur- 
rent channel is opened up and a jump in the I-V 
curve appears. The shape and height of these cur- 
rent steps depend on the width of the molecular 
resonances. With the increase of the molecule-to- 
clectrode coupling strength, the current varies con- 
tinuously with the applied bias voltage and gains 
much more bigger values, as shown in Fig. [5]^b). 
This behavior can be clearly observed by noting the 
area under the curve of the conductance spectrum 
given in Fig. ^h). Comparing the results plotted 
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Figure 5: Current / as a function of the applied bias 
voltage V for the phenalenyl molecule connected to 
the electrodes according to Fig. [IJa). (a) and (b) 
correspond to the results for the weak and strong 
molecule-to-clcctrode coupling limits, respectively. 

in Figs. [Sl^a) and (b) it is clearly observed that the 
current amplitude gets an order of magnitude en- 
hancement with the increase of the molecular cou- 
pling strength. So the electron transport through 
the molecular bridge is significantly affected by the 
molccule-to-clcctrode coupling strength. 

The effect of the molecule-to-electrode interface 
structure on the electron transport through the 
phenalenyl molecule is much more clearly explained 
from Fig. [6] where we plot the currents for the 
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molecular bridge given in Fig. [2ljb). Figures a) 
and (b) correspond to the currents in the two lim- 
iting cases as in Fig. [S] Similar to the previous 
case, here also the current amplitude gets an order 
of magnitude enhancement with the increase of the 
molecular coupling strength. But the significant ob- 



pear. For the case t/U « 1, the resonances and 
anti-resonances would split into two distinct nar- 
row bands separated by the on-site Coulomb energy. 
On the other hand, for the case where t/U >> 1, 
the resonances and anti-resonances would occur in 
pairs. 
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Figure 6: Current / as a function of the applied bias 
voltage V for the phenalenyl molecule connected to 
the electrodes according to Fig. [Hb). (a) and (b) 
correspond to the results for the weak and strong 
molecule-to-electrode coupling limits, respectively. 

servation is that for this bridge (Fig. [2l^b)) the cur- 
rent amplitudes get reduced enormously compared 
to the results predicted for the molecular bridge 
given in Fig. [21(a) (see the results given in Fig. [5]). 
This is solely due to the quantum interference ef- 
fect between the different pathways that the elec- 
tron can take. So it can be emphasized that design- 
ing a molecular device is strongly influenced by the 
molecule-to-electrode interface structure. 

All the above mentioned characteristics are also 
valid if the electron-electron interaction is taken 
into account. The main effect of the electron corre- 
lation is to shift and to split the resonant positions. 
This is due to the fact that the on-site Coulomb re- 
pulsive energy U gives a renormalization of the site 
energies. Depending on the strength of the nearest- 
neighbor hopping integral {t) compared to the on- 
site Coulomb interaction ([/) different regimes ap- 



4 Concluding remarks 

To summarize, a parametric approach based on 
the tight-binding model has been used to inves- 
tigate the electronic transport characteristics of 
a single phenalenyl molecule sandwiched between 
two metallic electrodes. Here we have used the 
Green's function method to determine the trans- 
mission probability (T), the conductance {g) and 
the current-voltage {I-V) characteristics through 
the molecule. This technique can be used to study 
the electronic transport in any complicated molec- 
ular bridge system. This is an interesting molecular 
system to study because of the unique behavior of 
the phenalenyl molecule. In this article, contact- 
ing the phenalenyl molecule from two different lo- 
cations leads to a large change in magnitude and 
shape in the transmission spectrum and the as- 
sociated current-voltage curves. Electron conduc- 
tion through the molecule is strongly influenced by 

(a) the molecule-to-electrode coupling strength and 

(b) the interface structure. These findings indicate 
that designing a whole system that includes not 
only molecule but also molccule-to-electrode cou- 
pling and interface structure are highly important 
in fabricating molecular electronic devices. 

More studies are expected to take the Schot- 
tky effect, comes from the charge transfer across 
the metal-molecule interfaces, the static Stark ef- 
fect, which is taken into account for the modifi- 
cation of the electronic structure of the molecular 
bridge due to the applied bias voltage (essential es- 
pecially for higher voltages). Here we have also 
ignored the effects of inelastic scattering processes 
and electron-electron correlation to characterize the 
electron transport through such molecular bridge 
systems. 
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